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ABSTRACT 

We investigate the importance of ambipolar diffusion and Hall currents 
for high- resolution comprehensive ('realistic') photospheric simulations. To 
do so we extended the radiative magnetohydrodynamics code MURaM to 
use the generalized Ohm's law under the assumption of local thermodynamic 
equilibrium. We present test cases comparing analytical solutions with numerical 
simulations for validation of the code. Furthermore, we carried out a number 
of numerical experiments to investigate the impact of these neutral-ion effects 
in the photosphere. We find that, at the spatial resolutions currently used 
(5-20 km per grid point), the Hall currents and ambipolar diffusion begin to 
become significant - with flows of 100 m/s in sunspot light bridges, and changes 
of a few percent in the thermodynamic structure of quiet-Sun magnetic features. 
The magnitude of the effects is expected to increase rapidly as smaller-scale 
variations are resolved by the simulations. 



Introduction 



In recent years, three-dimensional numerical magnetohydrodynamics (MHD) 
simulations have been a boon to our understanding of a variety of phenomena in the solar 
atmosphere and underlying convection zone. By including an increasingly extensive array of 
important physical effects e.g. radiative transfer, magneto-convection, heat conduction along 
magnetic field lines, time-dependent hydrogen ionization etc such models are beginning to 
unravel many mysteries of solar magnetism. To varying degrees, the most im portant physi cal 



effects have been captured by state-of-the-art radia t ive MHD sinaulations (jVogler et al. 



2005 



Stein &: Nordlund 



2006 



Hansteen et al. 



2006 



Steiner et al 



2008 



Gudiksen et al 
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201ll ). which aim to model the solar atmosphere as realistically as computationally 



feasible. The ability of these models to reproduc e the observationa 
wide variety of p. 



lenomena, such as solar faculae ( 



De Pontieu et al. 



20071. 



20081: 



Stein et al.i 



20061). emerging flux regions ( Stein fc Nordlund 



Martinez- Svkora et al 



i2010 



20101), s unspots (IRempe 



Cheung et al. 



spicules (IHansteen et al 



et al 



2006 



2008 



2009a 



20101) . p ores (iCameron et a 



2009a b 



Keller et al 



2004' 



properties for a 



Car 



2006 



sson et al. 



2004; 



Cheung et al. 



Tortosa-Andreu &: Moreno-Insertis 



Cheung et a 



Martfnez-Sykora et al. 



2010 



2009b 



2007 



2009 



Kitiashvili et al. 



Stein et al 



201 ih and 



20111 ) etc is encouraging. 



However it is also clear that the magnetic diffusivities and viscosities used in the 
simulations are too high: the molecular values correspond to length scales much smaller 
than the presently achievable resolution . In these circumstances t he simulations have bee n 



moving towards ever 



2009 



Moll et al. 



20111 ). 



ligher resolution ( IVogler fc Schiissler 



2007 



Pietarila Graham et al. 



Indeed the simulations have reached the point where, at some locations in the 
atmosphere, we are approaching some of the lengthscales associated with the microphysical 
processes. In particular both the Hall current and ambipolar diffusion begin to influence 
the plasma dynamics in the photosphere at resolutions of 5-20 km and on timescales of 
seconds to minutes. The Hall effect results from the relative immobility of ions (relative 
to electrons) to react sufficiently quickly to the Lorentz force. Associated with this Hall 
drift of ions with respect to electrons is an electric field parallel to the Lorentz force. This 
effect is important for dynamical changes in the plasma at frequencies between the ion 
and electron gyrofrequencies {uj = eB/rrijc). For a fully ionized plasma and a magnetic 
field in the kiloGauss range, the proton gyrofrequency is Up ~ 10^ Hz and the Hall effect 
thus can be safely neglected. The solar photosphere, however, is only weakly ionized and 
the gyrofrequency quoted above does not apply. In a weakly ionized plasma, most of the 
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inertia is in neutrals, which do not respond directly to the Lorentz force. The neutrals 
do however feel the magnetic field indirectly through collisions with the ions. Since the 
photosphere is collisionally dominated, the neutrals effectively give the ions extra inertia, 
and the Hall effect becomes important at much lower frequencies. The relevant frequency 
becomes u* = eB/m*c, where m* = {gi + Qn)/ni is the effective ion mass. For the low 
values of ionization fraction in the photosphere and lower chromo sphere, the effective 



proton gyrofrequeii cy (or Hall frequency) can be as small as 10 Hz (jPandey fc Wardk 



Pandey et al 



20081) 



2008 



In addition to the Hall effect, ambipolar diffusion can have an important influence 
on the dynamics of the photosphere and chromosphere. Since neutrals do not directly 
experience the Lorentz force, they have the ability to drift across magnetic field lines. Like 
classical Ohmic diffusion associated with a scalar magnetic resistivity, ambipolar diffusion is 
accompanied by the dissipation of an electric current (the Pedersen current), which converts 
magnetic energy into heat. Unlike Ohmic (associated with finite Spitzer conductivity) and 
ambipolar diffusion , the Hall effect is not diss ipative and hence does not directly contribute 



to plasma heating (jPriest 



1982 



Parker 



20071 ). Nevertheless, the Hall effect does redirect 



magnetic energy throughout the system (i.e., it contributes a term to the Poynting flux) and 
thus potentially has important indirect influences on wave damping and plasma heating. 

Although state-of-the-art 3D radiative MHD simulations have so far neglected the Hall 
effect and ambipolar diffusion, there is a considerable body of theoretical and observational 
work that point to the importance of neutral-ion interactions in the solar atmosphere. For 



heatine have been investigated bv De Pontieu & Haerende 


1 ( 


1998). 


Goodman ( 


2000. 


2004. 


2005): Goodman & E 


[^azeminezhad 


( 


201C 




Khodachenko et al. 


(2004 




Leake et al. 


(2005 


), 


Pandev et al. 


( 


2008 




Sineh & Krishan 


c 


2010 


) and \ 


<!homenko & CoUadof 


5 (2011) 
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All these studies conclude that neutral-ion collisions (i.e. Pedersen contribution to Joule 
heating) is the dominant dissipation mechanism in the chromosphere (as opposed to a 
single-fluid MHD model with only Ohmic dissipation from Spitzer resistivity) . Furthermore 



the ensuing chromospheric heating may be the driver of spicules (jPe Pontieu fc Haerendel 



1998 



James et al. 



2003 



Erdelyi &: James 



20041 ). Sakai and co-authors have used two- fluid 



(neutral hydrogen and proton components) MHD simulation s to model the collision of 



current-ca rrying loops in t 



microjets (ISakai fc Smith 



le ch romosphere (ISakai et al 



20061 ) as well as to model penumbral 



20081 ). In the former study, they found substantial differences in 



the results between single-fluid MHD and corresponding two-fluid MHD simulations, from 
which they conclude that the single-fluid MHD approximation (for fully-ionized plasmas) is 
inappropriate for the description of chromospheric dynamics. Additionally, they note that 
their single-fluid MHD simulation neglects Pedersen dissipation (i.e. ambipolar diffusion), 
which should be included to capture the relevant neutral- ion effects. 



Leake fc Arberl (120061 ) and 



Arber et al. 



( 120071 ) performed idealized numerical 



simulations of magnetic flux emergence into the solar atmosphere. Using a Id height- 



dep end ent profile for t 



le hyd rogen ionization fraction consistent with the VAL-C quiet-Sun 



model ( Vernazza et al. 



1981 



y) they performed flux emergence simulations with and without 
ambipolar diffusion. A comparison between simulations with and without ambipolar 
diffusion showed marked differences. First of all, ambipolar diffusion preferentially 
dissipates the Pedersen current. Since the Pedersen current is perpendicular to the magnetic 
field, its dissipation leaves the field reaching the corona in a force-free state. Furthermore, 
the departure from the frozen-in condition of the magnetic field allows it to traverse the 
chromospheric layers with substantially less mass uplift (i.e. less work against gravity) than 
in simulations without ambipolar diffusion. 



While most of these theoretical studies have concentrated on the chromosphere, the 
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effects are also important in the upper photosphere. In hght of the potential importance of 
the Hall effect and ambipolar diffusion on solar atmospheric dynamics, the purpose of the 
present study is to investigate their influence on a number of physical scenarios in the solar 
photosphere. The paper is structured as follows: Section [2] describes model equations and 
implementation of the numerical code. Section [3] presents results from tests cases used for 
validation of the code as well as in obtaining physical insight into the basic properties of 
the effects. We then study the effect of ambipolar diffusion on the structure of a quiet-Sun 
magnetic flux sheet and the effect of Hall currents in a sunspot light bridge. Our conclusion 
is that we have indeed reached the point where their effects are becoming dynamically 
important. 



2. Equations and implementation 



The radiation MHD code MURaM flVoder et al. 



2005 



Rempel et al 



2009ah aheady 



takes into account important physical effects such as radiative energy exchange in the energy 
equation as well as partial ionization in the equation of state. Both are important for the 
energy balance and stratification of the convection zone and overlying photosphere. Here we 
extend the model to incorporate the Hall effect and ambipolar diffusion. In comparison to 
the MHD equations for ful ly-ionized pl a smas, this extension inv olves two additional terms 



in generalized Ohm's law (IParker 
equation is 



2007 



Pandey fc Wardle 



20081 ). The modified induction 



— = Wx[vxB-—lj jxB + ijxB)xB ], (1) 

dt C erie CQiVin 



where B is the magnetic field, Qi and Qn the ion and neutral mass densities, respectively, 
D = Qn/{Qi + Qn) is the neutral mass fraction, v = Dvn + (1 — D)vi is the bulk velocity of 
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the center of mass of the neutral and ionized components, j = c(47r)~^V x 5 is the current 
density, rj is the Ohmic diffusivity, c the speed of hght, e the electron charge, the electron 
number density, and the rate of ion-neutral collisions. The third and fourth terms on 
the right-hand side represent Hall effect and ambipolar diffusion, respectively. 

Eq. ([T]) can be written in the form 

dB ^ ^ 47r 

— = Vx[uxB-—r]j], (2) 

where 

U = + WHall + ^^Amb, (3) 

iTHaii = -HVxB, (4) 
VAmh = M(V xB)xB, (5) 

H = (6) 

Airene 

M = . (7) 

Here, u is the effective velocity operating on B. In the presence of the Hall effect and 
ambipolar diffusion, departures of u from v means that magnetic field lines are no longer 
frozen into the plasma. This consequence, however, is not equivalent to a change in the field 
topology. Since the two effects act like velocities advecting, compressing and stretching the 
field, the field topology is in fact preserved. Another way to look at this is to consider the 
electric fields associated with the two effects. Both the Hall and ambipolar electric fields are 
perpendicular to B, which means that they do not lead to magnetic reconnection. Ohmic 
diffusion must act in order for reconnection to take place. 

The MURaM code solves the MHD equations in conservation form, which means that 
physical quantities on the grid are updated based on the net fluxes of the quantities into 
and out of individual grid cells. This implementation ensures that the total mass, energy. 
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and momentum within the simulation domain are conserved in the absence of net fluxes 
through the domain boundaries. For studying how the Hall and ambipolar terms change 
the energy budget in the solar atmosphere, we implemented the effects in such a way such 
that the numerical scheme remains conservative. 

By virtue of Eqs. and ([3]), the generalized induction equation can be written in 

conservation form as follows: 

dB ^ -> 

+ V -[u® B - B ®u] = -V X {riV X B), (8) 

where u®B and B®u denote tensor products. The energy equation must also be modified. 
In conservation form, it reads 



de 
di 



+ V- 



-gv"^ + ge + p) + u 



B^_]_ 
Air Air 



B{u-B) 



V ■ {B X r]V X B) + V ■ (v-t) 



An 

+V ■ (KVT) + g{g ■ v) + Q 



rad; 



(9) 



where e = ge + \gv'^ + B'^ j'^ix is the total energy density, r the viscous tensor, e is the 
specific internal energy density, g the gravitational acceleration {g = 2.1 A x 10^ cm s^^ 
at the solar surface), K the thermal conductivity, and Qrad the radiative heating/cooling 
term, which is dete rmined by means o f solving the radiative transfer equations along rays in 



various directions (IVogler et al 



20051). The gas pressure p and temperature T are functions 
of g and e and are determined from a look-up table during runtime. In the absence of the 
Hall and ambipolar effects (i.e. u = v), Eq. ([9]) reduces to the energy equation solved 
in the original version of MURaM. The momentum equation remains unchanged since 
terms that are quadratic in the drift velocity vd = Vi — Vn (i.e. the difference in the 



bulk velocities of the ion and neutra' 



conditions (IHasan fc Schiiessler 



1985 



species) can be sa : 



Pandey fc Wardle 



'ely ri eglected under photospheric 



20081 ). 
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MURaM uses 4-th order Runge-Kutta explicit time integration. The inclusion of 
the Hall and ambipolar terms in the induc tion and energy equations iiitrodu ces new 



constraints on the CFL stability condition ( ICourant. Friedrichs. &: Lewy 



19281). Since the 



Hall and ambipolar effects express themselves in the equations as effective velocities v, the 
associated time-step constraints imposed by the CFL condition are AtHaii = Aa;/|'yHaii| and 
A^Amb = '^x/\vAmh\ where Ax is the grid-spacing. Since WHaii and VAmh are first-order in 
j ~ AB/Ax, both AtHaii and AtAmb are in fact quadratic in Ax. For this reason, both the 
Hall and ambipolar effects can be considered as non-linear diffusion terms in the induction 
equation. 



2.1. Coefficients in the Hall and ambipolar terms 



The coefficients H = c(47rene)~^ and M = D^/cQiUin in Eq. and ([7]) are 
precalculated onto a look-up table for use during runtime. There are several assumptions 
made in creating the tables. The first is that the plasma is in local thermal equilibrium. The 
second assumption is the relative fraction of each element in the plasma. The Saha-Eggert 
equations can then be solved to obtain the number densities of eac h chemical specie s and 
the number density of free electrons. We use the code described in IWittmannI (jl974j ). The 



coefficient in Eq. ([6]) can be determined without further assumptions. The coefficient in Eq. 



m 



requires an addi t ional approximation to obtain z/j„. Three different sources are given 



De Pontieu et al. 



Osterbrock 



fll96lh . 



(120011 ): we chose to use choice 'a' which corresponds to that given by 



In the chromosphere, the plasma is sufficiently tenuous such that U 



a suitable assumption for determining ionization fractions of hydrogen (IKneer 



E is no lon ger 



1980|). 



Therefore, a time-dependent app roach to calculating n umber densities of the different 



hydrogen specifies must be used (ILeenaarts et al 



20071 ). This is not yet implemented in 



the code. However, the numerical experiments presented below focus on dynamics at the 
photospheric layers. 



3. Test cases 

Here we describe two simulation setups used as test cases for validation of our 
implementation of the Hall and ambipolar diffusion effects. 

3.1. Influence of the Hall effect on plane-polarized Alfven waves 

Consider the one-dimensional problem of Alfven waves propagating along the 
x-direction in the presence of a guide field, B^. The transverse {By and B^) components 
of the field are assumed to be small compared to the guide field B^, so that second-order 
terms in these quantities can be neglected. For a pure Alfven wave, the total background 
mass density (^o) is constant and uniform. For this test example we neglect the effect of 
ambipolar diffusion. The evolution equations are then 



= iri^ xB)xB, and 

ot 47r 

OB r -» -> 

— = Vx vxB - H{W X B)xB 



(10) 



(11) 



Let the solutions for v and 5 be of the form 








y 




v'e 



;(a;t-M[cos at -id sin at] 
i^t-i^^)[smat + id cos at] j 
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y ^/gi(a;t-fex)sinat ] 



where a and 5 are still to be determined, v' and 6' are amplitudes for the traveling wave 
solutions and are assumed constant (for a chosen k\ The y-component of the momentum 
equation reduces to 



f '[iw cos at — a sin at + 5uj sin at — i^o" cos at] 



-ikbo 

-0 cos at, 



4:11 go 



(12) 



Since f' and b' are constants, equating coefficients of sin at in Eq. (fT2!) gives a = ^w. 
Substituting this back into Eq. fll2p and eliminating the common factor cos at, we obtain 



u kbo 

b 



iJ2 _ ^2 47]-^Q 



(13) 



Now the y-component of the induction equation gives 



b'[iuj cos at — a sin at] = ~ikbQv'[cos at — iS sin at] — b'Hbok'^ sin at 



Using Eq. (fT2!) to substitute for v' in ([H]), we obtain 



(14) 



iu cos at — a sin at = i 



blk^ 



[cos at — Z(5 sin at] — Hb^k'^ sin at (15) 



Equating coefficients for cos at in Eq. (fT5l) . we obtain the dispersion relation 



2 2 ; 2 I 2 

u = Vjsk + a , 



(16) 



where va = &o(4vrf)o) is the Alfven speed. Similarly, equating coefficients for sin at in 
Eq. f|T5l) and using the dispersion relation flT6|) gives 



a 



Hbok^ cbo 
2 STrrZee 



(17) 
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The dispersion relation Eq. f lT6|l is similar to that of a pure Alfven wave except for a 
correction term. For the phase speed u/k, the correction term is of order Od-UHaiil)? where 
t'Haii is the Hall velocity as defined in Eq. (jl]). The physical interpretation of the above 
result is as follows. Suppose the Alfven wave of wavenumber k is initially plane-polarized 
in the y-direction. Under the action of the Hall effect, the plane of polarization of this wave 
would precess about the at a rate equal to a. 

For the following one- dimensional test case, we choose to use physical parameters that 
are typical of photospheric conditions. The a;-component of the magnetic field is uniform 
with B.j.{t = 0) = 60 = 100 G. Superposed on this uniform field is a sinusoidal profile of 
the tranverse field By{x,t = 0) = b'cos{kx), where b' = 0.1 G and k = 2tt x 10^^ km~^. 
The one-dimensional simulation domain has width = 100 km so the domain is spanned 
by exactly one wavelength of the transverse perturbation in By. The initial total pressure 
of the unperturbed magnetic field and plasma is ptot = Pgas + -B^/Svr = 4 x 10^ dyn cm~^. 
The mass density is approximately g = 10^'' g cm~^ and the temperature is T = 6000 K. 
Gravity is switched off so that the only force initially acting on the plasma is magnetic 
tension. For the sake of carrying a test simulation, we impose a constant Hall coefficient of 
H = 10* cm^ s^^ G^^ and switch off the ambipolar term (i.e. M = 0). The grid spacing 
used for the simulation is Ax = 1 km and periodic boundary conditions are used. 

The black and purple dots in Fig. [H respectively, show values of By and Bz at x = 
sampled at different times during the simulation. The black and purple curves in the same 
plot show the corresponding analytical solution to the problem. Since the initial condition 
consists of the standing wave, the solution consists of a superposition of modes propagating 
in the positive and negative x-directions, viz. 

b' 

By{x,t) = — cos{at)[cos{kx — cut) 

— cos{—kx — ut)], (18) 
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Time [sec] 



Fig. 1. — The two magnetic field components By (black) and (purple) as functions of 
time (sampled at x = 0) from the numerical (dots) and analytic (solid line) solutions for 
k = 27C X 10~^ km~^. The fast time variation corresponds to the Alfven waves, the slow 
variation to a rotation of the plane of polarization due to the Hall effect. 




Tirpie [sec] 



Fig. 2. — Same as Fig. [T]but with twice the wavenumber. Note that the time interval shown 
in this plot is shorter than that shown in Fig. [H 
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= b' cos{at) cos{kx) cos{ijjt) , (19) 
Bz{x,t) = b' sm{at) cos{kx) cos{ut) . (20) 

For the set of physical parameters given above, the rate of precession of the plane 
of polarization is cr = 2 x 10~^ Hz. This corresponds to a rotation period of 
r = 27ia^^ = 3.2 x 10^ s. The intrinsic oscillation frequency of the Alfven wave is 
u = 5.5 X 10^^ Hz. For the same initial setup, but with twice the wavenumber (i.e. 
k = Arc X 10"^ km), oo = 0.11 Hz and cr = 8 x 10~^ Hz. Fig. [2] shows the solution for 
this case. Compared to the case when k = 27t x 10^^ km^^ (Fig. [II), we see here that for 
k = Att X 10"^ km^^ there are half as many Alfvenic oscillations for each complete rotation 
of the plane of polarization since a oc k"^ and u oc k. 



3.2. Collapse of current layers due to ambipolar diffusion 



An effect of ambipolar diffusion is the reduction of the width of current sheets (jParker 



1963 



Brandenburg fc Zweibe 



1994J ). In the absence of Ohmic diffusion, the time- 
independent solution for the transverse magnetic field in the neighborhood of the tangential 
discontinuity has the profile 

B(xx^/\ (21) 

where x is the coordinate in the direction normal to the curernt sheet. They also point out 
that the current density at the neutral line increases without bound if Ohmic diffusion is 
absent (namely, the constant of proportionality in Eq. fl2T]) increases with time). Acting 
alone, ambipolar diffusion does not lead to magnetic reconnection. This property is a direct 
consequence of the fact that the ambipolar term in the induction iquation can be written as 
an effective velocity "UAmb acting on B. So ambipolar diffusion merely acts to sharpen the 
current layers which in turn allows Ohmic diffusion to occur at an enhanced rate. 
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Fig. 3. — The magnetic field (upper panel) and current density (lower panel) distributions 
for the test case corresponding to a collapsed current sheet in which Ohmic and ambipolar 
diffusion balanced each other. In both panels, the solid line indicates the semi-analytical 
solution and the purple dots indicate the solution from the simulation. The deviation of the 
two away from the current sheet is expected. 
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To test our implementation of ambipolar diffusion, we consider a one-dimensional 
system with the following initial distribution: 

B{x,t = 0) = {0,Bosm{kx),0). (22) 

We set Bq = 1000 G, k = 2tc x 10~^ km~^ and x G [0, 100] km. The plasma has an initial 
uniform temperature of T = 6000 K and the system is in total pressure equilibrium, with 
Ptot = Pgas + B'^/SiT = 4.4 X 10^ dyn cm~^. For this simulation, we set a constant Ohmic 
diffusivity of ?7 = 10^ cm^ s"^, a constant ambipolar coefficient of M = 5 x 10^ cm^ s~^ 
G^^, and switch off the Hall term (i.e. H = 0). Periodic boundary conditions at x = and 
X = 100 km were imposed. 



We follow the analysis of iParkerl (Il963[ ) to seek a time-independent solution for the 
collapsed current sheet including Ohmic diffusion. For the one-dimensional setup described 
above, the induction equation reduces to 



dBy _ d 
dt dx 



{MBl + rtf^^ 



(23) 



dx 

For a time-independent solution and for constant M and rj, Eq. ( I23i) can be integrated to 
yield 



M 



g Bl + r^By = Cx + D, (24) 

where C and D are constants determined by boundary conditions. For our setup here. By 
is antisymmetric about a; = 0. This means By{x = 0) = 0, which implies D = 0. Since the 
simulation setup employs periodic boundary conditions at a; = and at x = 100 km, there 
is not an infinite reservoir of magnetic fiux which can be transported into the diffusion 
region for reconnection. As such, the simulted system never really reaches a true steady 
state and the integration constant C depends on how much flux has reconnected from the 
initial state. 

In the absence of Ohmic diffusion [t] = 0), Eq. can be directly solved to yield 
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the relation given by Eq. fl2T]) . Otherwise, we identify that for a specific spatial location x 
(distance from the neutral line), Eq. fl2^ is a cubic polynomial in By. To determine for 
By{x) we used a numerical solver for the polynomial corresponding to various values of x. 
This is plotted as solid lines in both panels of Fig. 13.21 The corresponding solution from 
the simulation is plotted as purple dots. For C = 2.85 x 10^ G cm s^^, the simulated and 
semi-analytical solutions closely resemble each other in the vicinity of the current sheet. 
For reasons already discussed, they are expected to deviate far away from the current sheet. 

4. Simulation results 

In the following sections we present simulation results on how the additional effects 
impact the physical evolution of magnetic field in conditions mimicking the solar atmosphere. 

4.1. Magnetic reconnection in the presence of ambipolar diffusion and 

radiative cooling 

In Section [3?2| we examined how ambipolar diffusion steepens gradients in the magnetic 
field to sustain current sheets. For the sake of code validation, we assumed a constant 
ambipolar coefficient M. However, this assumption neglects an important negative feedback 
loop, which proceeds in the following fashion. Both ambipolar diffusion and Ohmic 
dissipation act to heat plasma in and around the non-ideal region. The heating of the 
plasma leads to an increase in the ionization degree of the plasma, which in turn could 
suppress the efficiency of ambipolar diffusion (i.e. decrease M). So although ambipolar 
diffusion may act to steepen up a current sheet, it does not necessarily result in runaway 
magnetic reconnection. 

To investigate this effect, we carried out three one-dimensional numerical experiments 
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(see Table [T]). The initial condition in all cases consists of a current sheet with 

By = 1000 G, if X > 0, 

= -1000 G otherwise. (25) 

The other components of B are zero. The plasma initially has uniform temperature of 
T = 5000 K and the uniform total pressure ptot = Pgas + B'^/Stc = 4.13 x 10^ dyn cm~^ 
(except at the tangential discontinuity). This choice leads to (3 = Sirpgas/ B"^ = 0.04, so that 
dynamics are dominated by magnetic forces. The domain spans x G [—300, 300] km and the 
grid-spacing is Ax = 6 km. Periodic boundary conditions are imposed at x = ±300 km. 

For the initial thermodynamic state of the plasma, Eq. ([7]) gives M = 2.3 x 10^ cm^ 
s-i G-\ In run ReA, we keep M uniform and constant with this value. In run ReB, we 
allow M to have values consistent with the local thermodynamic state at any time by using 
our look-up table. The radiative cooling term in the energy equation is switched off in 
both runs ReA and ReB. In run ReC, we introduce a cooling term in the energy equation 
to mimic the energy loss from the plasma as it heats up above the temperture of its 
surroudings. MURaM is capable of calculating this source/sink term self-consistently by 
solving the radiative transfer equation along multiple directions, but this is not appropriate 
for a one- dimensional setup. The cooling term in run ReC is of the form 



Run 


Ambipolar coeff. 


M 


Cooling 


Flux cancelled by t = 820 s 


ReA 


2.3 X 10^ cm^ s-i 


G-2 


No 


28% 


ReB 


Variable 




No 


24% 


ReC 


Variable 




Yes 


51% 



Table 1: Set of reconnection experiments with and without variable ambipolar diffision and/or 
Newton cooling. Eq. ([7]) is used for evaluation of the variable ambipolar diffusion M. 
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1 .0 




-200 -100 100 200 
X [km] 



Fig. 4. — The distribution of magnetic field, temperature and ambipolar diffusivity r/Amb in 
simulation runs ReA (solid line), ReB (dashed line) and ReC (diamonds) at time t = 820 s. 
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Qr.d = -ag[T^-T^], (26) 

where g is the local mass density, T the local plasma temperature, Tq = 5000 K and 
a = 10~^ erg s~^ g"^ K""^. T > Tq implies Q^ad < and the plasma loses energy. This 
decreases the plasma temperature (and ionization fraction) and provides a way for the 
ambipolar diffusion to continue acting. 

Figure m plots the distribution of By, T and the ambipolar diffusivity r^Amb = MB"^ for 
all three simulation runs at time t = 820 s. Between t = and this instant, the fraction of 
the initial unsigned flux that has cancelled in runs ReA, ReB and ReC is 28%, 24% and 51%, 
respectively. The reason more flux has cancelled in run ReA than in run ReB (a difference 
which increases with time) is that in the latter, the temperature in the current sheet due 
to Ohmic and ambipolar dissipation suppresses M and thus ?7Amb- This negative feedback 
mechanism is absent in run ReA since the ambipolar coefficient is set to be constant. In 
run ReC, radiative cooling in the current sheet maintains the action of ambipolar diffusion 
and leads to a sustained flux cancellation rate of almost double that of the two other runs. 

4.2. Vertical flux sheet 

We have also investigated the effect of ambipolar diffusion on the structure of 
small-scale magnetic features in the photosphere. For this purpose we set up a two 
dimensional simulation extending vertically from about 700 km below to a height of 700 
km above the photospheric base (tross = !)• The horizontal extent of the box was chosen 
to be 500 km and contains one convective roll. This was done so that only one flux sheet 
is ever present which simplifies the analysis. A grid spacing of 5 km was used in the 
horizontal direction and 14 km in the vertical direction. Non-magnetic convection was 
allowed to develop for several hours, after which an initially uniform vertical field of 100 G 
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was imposed, and the system was evolved for several additional hours. This state was then 
used as an initial condition for runs with and without ambipolar diffusion. Fig. \5\ shows the 
temperature, density and vertical component of the magnetic field after 2.5 minutes for the 
case without ambipolar diffusion and the difference resulting from the inclusion of ambipolar 
diffusion. The effect is relatively weak (a few percent) consistent with expectations (e.g. 



Pandey fc Wardle 



20081 ). 



4.3. Magnetoconvection in a sunspot umbra 



In this simulation, we examine the influence of the Hall effect on magnetoconvection 
at field strengths found within sunspot umbrae and in the vicinity of sunspot light bridges. 
In this two-dimensional simulation, we begin with statistically relaxed hydrodynamical 
convection driven by radiative cooling at the photospheric base. The computational domain 
spans 18.4 Mm in the horizontal direction and 6.1 Mm in the vertical direction with 
horizontal and vertical grid spacings of 9 km and 6 km respectively. After introducing a 
uniform vertical magnetic field with = 1500 G, we ran the simulation for approximately 
one hour to allow the magnetoconvection to relax. The result is the development of 
strong magnetic regions with deep Wilson depressions interspersed with upfiow regions 



with relatively weak fielc 



dots (iSchiissler fc Vogler 



This type of cross-sectional st ructure may represe nt umbral 



20061 ) and sunspot light bridges (ICheung et al.l 120101 ) 



Figure |6] shows the distribution of physical quantities of a structure in the simulation 
resembling a sunspot light bridge. In both the temperature and vertical magnetic field 
panels, one fin ds the characteristic cusp shape where canopy fields from adjacent umbral 



regions merge (ISchiissler fc Vogler 



20061 ). In the absence of the Hall effect, the magnetic 
and velocity fields are confined to the x- and 2;-components. The sharp interfaces between 
weak-field regions within the light bridge and the strong field regions outside it acts as the 
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Fig. 5. — Effect of ambipolar diffusion on the evolution of a small magnetic flux sheet in 
the photosphere. The left panels show snapshots of the vertical component of the magnetic 
field (top), the mass density (middle) and the gas temperature (bottom) for the run with 
ambipolar diffusion switched off. The right panels show the difference between the snapshot 
and a similar snapshot where ambipolar diffusion was included. These snapshots are taken 
after 2.5 minutes of evolution from a common initial condition. 
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Fig. 6. — Cross-sections of physical quantities from a 2D simulation of umbral magnetocon- 
vection with treatment of the Hall effect. The Hall term leads to an out-of-plane component 
of the magnetic field {By, lower left panel). The presence of this third component leads to a 
Lorentz force which acts to accelerate plasma in the out-of-plane direction {vy, lower right 
panel) . 
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site of high current density, which in the presence of the Hall term leads to an effective vnaii 
in the y-direction. This component of the effective velocity in the y-direction is non-uniform 
and leads to the generation of a y-component (out-of-plane) of the magnetic field. After By 
is generated by the Hall term, convective flows advect the flux from the layer and back into 
the convection zone. These flux rolls can be seen in the bottom left panel of Fig |6l 

We emphasize that the strength of By is only ~ 5 G, which is much weaker than the 
kilogauss fields in the umbra. Thus in the current simulation, the fields generated by the 
Hall currrent do not strongly influence the dynamical evolution of the system. Nevertheless, 
the presence of non-zero By leads to a Lorentz force which accelerates an actual (i.e. not 
just ■UHaii) out-of-plane flow Vy (see bottom right panel of Fig [6]) to speeds of up to 0.1 km/s. 
Such flows are completely absent in a simulation without the Hall effect. 



5. Discussion 

Both ambipolar diffusion and the Hall current have now been included in a code 
designed to produce 'realistic' simulations of the solar photosphere. The implementation 
has been validated, and the physical effects such as the rotation of the plane of polarization 
of Alfven waves and the steepening of current sheets has been reproduced. 

In simulations of magnetic reconnection with ambipolar diffusion, we find that 
ambipolar dif fusion acts to sharpen c u rrent layers. This is a well-known result from 



the study by [Brandenburg fc Zweibell ( 1l994j ). However, our study also shows that the 
temperature increase associated with ambipolar and Ohmic dissipation in the current sheet 
rapidly suppresses ambipolar diffusion. In this sense, ambipolar diffusion is a self-limiting 
mechanism for enhancing magnetic reconnection. Only when we include an energy loss 
mechanism for the plasma (in this case a radiative cooling term) does the effect remain 
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important. This result has important imphcations for studies of chromospheric dynamics. 
The effect of ambipolar diffusion in the chromosphere may be severely overestimated when 
the ionization fractions and energy loss mechanisms in chromospheric pla sma are not 
treat ed in a realistic fashion (e.g. such as time-dependent H ionization, see 



Leenaarts et al. 



20071 ). 



In a two-dimensional simulation of umbral magnetoconvection, we examined how the 
Hall effect generates a third (out-of-plane) component of the magnetic and velocity fields 
at interface layers between weakly magnetized light bridges and neighboring strong field 
umbral regions. At the grid spacing of ~ 10 km, the Hall term results in generation of the 
third component of the magnetic field at the level of ~ 5 G from kilogauss umbral fields. 
The generation of the out-of-plane component of the magnetic field faciliates a magnetic 
tension force which accelerates plasma out of the plane to speeds of up to 0.1 km/s. 

At the current grid resolution, the influence on dynamics of the system due to the 
presence of the Hall effect begins to become important. Furthermore, we expect that the 
relative importance of both the Hall and ambipolar effects to increase with grid resolution. 
This is due to the expectation that with decreasing grid spacing Ax, numerical simulations 
would be able to resolve increasing sharp current layers (since j ~ AB/Ax). Since both 
the Hall and ambipolar velocities (■yHaii and VAmh, respectively) in the induction equation 
are linear in j, their amplitudes relative to the actual fluid velocity v and to numerical 
diffusion should increase linearly with grid resolution. A Hall- induced velocity of 0.1 km 
s^^ at Ax = 20 km, when scaled to Ax = 2 km, already yields a velocity of 1 km s^^. Thus 
the hnear trend is unlikely to continue down to Ohmic dissipation scales (1 < Axohm ^ 10^ 
m in the photosphere). With this in mind, future investigations will examine how Hall 
amd ambipolar effects will, at sufficiently high spatial resolution, come to be important for 
dynamics at small scales. 
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